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ABSTRACT 


Eimenawake data from California from the period of 
ieamiaryel.) 1969, to December 31, 1971, has been analyzed 
byetwo methods. The first method compares the occurrence 
@emearthquakes with the Synodic, Draconic and Anomalistic 
MimtaeecrrOds..) otatistical tests to determine if the earth- 
G@uakes are uniformly distributed over each of the three 
Mitiarepereods are presented. The second method compares 
Eiree Components of Mi mcm TOrecavcetor and thei7 


Semivataves at the times of earthquake occurrences. 
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Poe NU RODUCT ION 


Piin@ecorre lation with the time of the occurrence 
of earthquakes as reported by Allen [Ref. 1] has been one 
@f the methods used in trying to determine the causes of 
weweengqiakes for the better part of a century. Knopoff 
Ref. 3] reported that it may be that earthquakes have a 
Peeater tendency to be triggered when the sun and the moon 
Baeenecarhy aligned with the earth. An interesting observa- 
tion which was reported by Wigand [Ref. 16] was that the 
Sun and the moon were almost aligned with the earth at 
miemtime Of the San Fernando Valley Earthquake of 
February 9, 1971. It is further contended by Wigand 
[Ref. 16] that there is strong evidence that some earth- 
feikes are related to tidal force, and some are actually 
imeueeered Dy the tidal force. The objective of this 
thesis was to investigate these findings by analyzing the 
Pmavitational forces of the moon and the sun with relation 
momrecorded carthquakes. 

Two approaches in analyzing earthquake data are 
eescmecd es ie fitnst approach was to compare the occurrence 
Peeeanmenquakes with the Synodic Lunar Period, the Draconic 
Ptimeereriod, and the Anomalistic Lunar Period. The second 
approach was to investigate the tidal force exerted by the 
sun and the moon with relation to 155 recorded earthquakes 
Whreh occurred in the San Fernando Valley from April, 1971, 


bomnevemocr, 1971, 





iitmeannecrlinar periods discussed in the first approach 

are as follows: 

Mo noute sLUnar Period, mean value 29.530589 solar 
days. Defined as new moon to new moon, where the 
HeCeOoWmmOoccunsSewhen the geocentric longitudes of 
the sun and the moon are the same [Ref. 12]. 

mee liaconic munar Period, mean value 27.212220 solar 
tipeweicrined ds tie ascending node to the 
ascending node, where the ascending node pas- 
TIemoaeene MOO, 1S from minus decrees latitude 
to positive degrees latitude [Ref. 12]. 

ae wonmalistic Lunar Period, mean value 27.554531 
Tmo sw DetLinecdad aS perigee to perigee, 

Witeme wer icce 15 that point of the moon’s orbit 

at which the moon is nearest to the earth [Ref. 12]. 

Pivmlunar pervods |Ret. 12) actually exist but the 

@pove three were included for the reasons given below. 
In analyzing the earthquake data with respect to the 
Palocie Linar Period, the occurrence of earthquakes with 
feespectneo Lie positions of the sun and the moon ae 
exauemea., = In analyzing the data with respect to the 
Draconic Lunar Period, the occurrence of earthquakes with 
neepeet to the moon's position in the Northern Hemisphere 
and Southern Hemisphere were examined. This may have been 
important since the data used included only earthquakes 
fiPmeeteoceurred in Calitornia. In analyzing the earthquake 


data with respect to the Anomalistic Lunar Period, the 





fee rence I) the gravitational force of the moon on the 
€arth was examined as the moon proceeded from perigee to 
peg@gee to perigee. 

itiemearchquake data received from the National Center 
ime@meeeartniquake Research spanned three time periods and 
Paetienhod 1m CwO areas Of California. The first set of data 
Peematedws lass Carthauakes recorded from April 1, 1971, to 
November 17, 1971. These are classified as aftershocks of 
the February 9, 1971, San Fernando Valley Earthquake 
maeoure 1)> The magnitudes ranged from 1.0 to 4.2 as 
Mmeaowred on the Richter Scale. This set of data will be 
meascnrred £0 aS the San Fernando Valley Data (SFVD). 

iiemsceond Set Of Gata included 2879 earthquakes 
meconced from January 1, 1969, to August 23, 1970. These 
earthquakes occurred in the San Francisco Bay Area (Figure 
Mein the vicinity of the San Andreas Fault, the Hayward 
Bowie, and the Calaveras Fault. The magnitudes ranged 
from 0.11 to 5.06 as measured on the Richter Scale [Refs. 
feomcdes|.  bhis set of data will be referred to as the San 
Francisco Bay Area Data I (SFBAD-I). 

iiecmeiindmseteor data included 2734 earthquakes recorded 
from August 1, 1970, to December 31, 1971. The area was 
the same as SFBAD-I. The magnitudes ranged from 0.29 to 
#.45 as measured on the Richter Scale [Refs. 5 and 6]. 
Paaisewset of data will be referred to as the San Francisco 


Bay Area Data II (SFBAD-II). 
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SFBAD-I and SFBAD-II both included earthquakes which 
occurred in August, 1970. In the analysis, as discussed 
in the next section, the duplicate earthquakes were eliminated. 

In analyzing the data with respect to the lunar periods, 
the following hypotheses were investigated: 

Null Hypothesis: Earthquakes occur in a 
Unaeeowinc! Serroue lon. wii Ch 
Reiee mE OmmtNe wlunay period . 
Alternative Hypothesis: Earthquakes do not 
COCCUre Iiea Uniform dilstri- 
DUtION with respect to the 
I ybeveie jorese alfexat - 

Two tests were used in the analysis: The Chi-Square 
Memietor olenificance [{Ref. 17], and the Goodness of Fit 
Statistic Vn [Ref. 13]. The tests are discussed in what 
follows: 

itiemena-square Test for Significance 

N, the total number of observations in the 
samples which were analyzed. 
H, the total number of equal intervals in 


which the lunar period was divided. 


X = (ae the expected number of observations which 
occurred in each interval of a uniform 
Ge Satie a Dulte Won: 
Ti, the total number of observations recorded within 


Cxlchmmiterval. 1=1,,..,H. 


10 





eile thne degrecs Of freedom associated with 


ile tee sin 
: H (T.-x)¢ 
y, Be > een mtn mOMl-oduare Statistic. 


jrl 


itiemGoodaness Of Pit Statistic Vn: 


N. the total number of observations in the 


sample analyzed. 

Tooele mUntironml distribution on the lunar 
period in which the origin was the origin 
Cumeiemluinan period , 

FyQd), Pema roerTONiTOn Ob tne ObSerVations on 
the lunar period with the same origin as EX) 

X, the time of each earthquake in the sample. 

ie oUpeCs COCO) ) = Int (F(x)-FOx)), 


-O<X<o@ -0O< xX <a 
, the 


Vieseatistic. 


Secpicicm | het, IS) stated that 1f the observations 
PeempOllts om a circle, the value of Vn obtained from the 
moovencaduation does not depend on the choice of the origin 
mouemeasurianeg x. <The distribution of Vn refers to the dis- 
Poenwatton under the null hypothesis [Ref, 13]. Tables of 
eomiticance Of the Statistic YN Vn are included in 
heierence 13. 

A simulation of 1000 earthquakes over a period of six 


months was included in the analysis. The simulation was 


EA 





ae oniledsuonbesu tne Uniform distribution of the events 
Sect une three lunar periods. 

The final analysis in this thesis was made possible 
iyeeenes use Of a Computer program by Professor Rex H. 
Shudde of the Naval Postgraduate School. This computer 
program calculated three components of the tidal force 
Meedsured if microgals) at the epicenter of an earthquake. 

The three components are the radial component, the 
north component, and the east component. The radial com- 
Poletti thestidal force vector is through the center of 
Micmcabuleandethce e€picenter of the earthquake. The north 
midmecast Components of the tidal force vector are in the 
milane which passes through the epicenter and is perpendicular 
to the radial component. The computer program also calcu- 
fated the rates of change of the three components described 
above. The three components were plotted for a 24-hour 
time period on the day each earthquake occurred. A visual 
Somparisom of the plots of the 155 earthquakes which 
Pectuueds im the Sdn fernando Valley with respect to the 
three components of the tidal force was made. The results 


CeeenC COlparison are im Section III. 


le 





MeeeenneLyols AND DISCUSSION OF THE LUNAR 
PERO pe COM? UTR PROGRAM 


The computer program used to analyze the frequency of 
earthquake occurrences during the Synodic, Draconic and 
Anomalistic Lunar Periods contains the following steps. 
iiemitrst Step involves the time frame of the three sets 
of data and the lengths of the lunar periods within these 
tame fLrames. the mean times of the three lunar periods, 
as given in the Introduction, were not used in the progran, 
Since the actual lengths of the periods vary considerably. 
An example of this difference was the Anomalistic Month 
Pimeiieiwas a mean period of 27.32166 days. During 1971 the 
length of the Anomalistic Month ranged from 25.41667 
Gays) tO 28.45855 days [Ref. 9]. 

The times of new moon and lunar perigee for the years 
moos L970 and 197] are tabulated in References 7, 8 and 
9. From these times (to be referred to as origins) the 
diafterent period lengths of the Synodic and Anomalistic 

Wit wom vere computed, —keterences 7, & and 9 also include 
Eaokeseindicating the apparent latitude of the moon for 
Miem—nreSpeetive years at intervals of one half day. The 
origins for the Draconic Month, as the moon passed from the 
southern Hemisphere to the Northern Hemisphere, were 
computed from these tables by linear interpolation. This 
peocceureesshould be fairly accurate since it involves only 


feces ana seconds of latitude. 


13 





The origins of the three lunar periods and the times of 
each earthquake occurrence were put into the program in 
year, month, day, hours, minutes and seconds in Greenwich 
Mean Time (GMT). Within the program the year, month and 
Mey Of Gach event were converted to the respective Julian 
Day Number (JD) [Ref. 14]. The hours, minutes and seconds 
Gemstituted a fraction (Fract) of a eye The Julian Date 


was computed as follows: 
wwtan lates JN + Fract - 0.50 


The one half day was subtracted due to the differences 
in origins of the Julian Day Number and GMT. The Julian 
Day Number had its origin at noon whereas GMT has its origin 
meeinidnight. 

In order to simplify the statistical tests which followed 
in the program, the data was truncated to insure that no 
fractions of a period were included. An example of this 
was the SFBAD-J and the Synodic Month. The time the 
SFBAD-I started was January 1, 1969, which was between the 
two new moons of December 19, 1968, and January 18, 1969. 
Since the events that may have occurred from December 19, 
foODoemtouvecemmper SF. 1908, were not included, the events 
from January 1, 1969, to January 18, 1969, were excluded 
in the analysis. Similarly fractions of periods were 
eliminated for the Draconic and Anomalistic Months. 

Each of the three lunar time periods was reduced to the 


unit circle as explained in the following example. The 
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temeeieg: the oynodic Lunar Period from April 25, 1971, 
momiay 24, 1971, was 29.354170 days. At the time of the new 
HeomeOnmapct! 25, the yalue on the unit circle was zero. 
Pemene time Of the first quarter moon on May 2, the value 
onmeme Unit Circle was 0.25. At the time of the full moon 
open) 9, the value on the unit circle was 0.50. At the 
time of the third quarter moon on May 17, the value on the 
Mamereircle was 0.75. Finally at the time of the new moon 
imi 24. the value on the wnit circle was 1.0 [Ref. 9]. 

The time at which each earthquake occurred during the 
itmar period (residual) was placed upon the unit circle 
mabe Gespecteto the respective origins. This was accom- 
plished by dividing the difference between the Julian Date 
of the earthquake and the Julian Date of the origin immedi- 
pwcly preceding the earthquake by the length of the period 
in which the earthquake occurred. 

Tie unit Circle was divided into Z0 equal intervals 
Mmondensrospcrtorm the Cha-square Test for Significance 
onl jeassexplained in the Introduction. Ostle [Ref. 10] 
Suggests, as a rule of thumb, to use no fewer than three for 
miemexpecreds frequency ameeach interval. Ostle [Ref. 10] 
Piernerec@ntends that if some Expected numbers are too small 
Picmena-square Statistic will be a poor indicator of the 
WimMclevOLe tne Mypornesis Under test, By selecting 20 
intervals, the earthquakes in the three sets of data comply 


MiiMemicmrUleceOn thumo £Or Expected frequencies. 
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litemrestaulals that were computed were then sorted into 
peecnaimesorder |Ref. 2} so that the Vn Goodness of Fit 
lestecould be computed as explained in the Introduction. 
The Vn Goodness of Fit Test is a modification of the 
ROlmogorov-Smirnov Test {[Ref. 13] which is considered to 
tmeemeraliy a more powerful test than the Chi-Square 
estm(| Ref. 10]. 

iijemcueut of the program consisted of the total number 
@mmpevencs that occurred during each of the three lunar 
Periods; the number of events and the percentages which 
meeurred within each of the 20 intervals, the Chi-Square 


statistic, and the Vn and YN Vn statistics. 


1:6 





Pee Rous 


MicmorMirauion or the 1000 earthquakes as stated in the 
introduction was performed to test the uniform distribution 
Oeene events on the three lunar periods. The null hypothe- 
Sis that the simulated earthquakes were distributed uniforn- 
iyewas accepted for each of the lunar periods. The bases 
fee the acceptances were the Chi-Square Test with 19 
Meprees Of freedom and the Vn Test. The Chi-Square statis- 


tics and the VN Vn statistics calculated are given below: 


Synodic Draconic Anomalistic 


mecha Souare Foes gaa 22.04 Ins Oe 
YN Vn iL, GON). e250 59 S:5 


Meenwiieiypothesas and ailltermative hypothesis are: 
MueienioOGheSiS se sarthquikes occur in a 
Di tOonmimarcotr lbUutton with respect 
Eoutnes linear sper Lod: 
Mesci Mae VewMMOGtneSTS scar encuakes do not occur 
in a uniform distribution with 
Pecpeetmuowume lunar period. 
ifeetesnuli hypothesis 1s accepted, at might imply that 
tt micComucmeravitational effects On the earth do not trigger 
earthquakes. Conversely if the null hypothesis is rejected 
and the alternative hypothesis is accepted then it might 
Mp weit the moon's eravitational effects on the earth 


iemeiregey Sarthaquakes. 


vay 





inewresults of the statistical tests performed on the 


Gitweewccts Of data wath respect to the Synodic, Draconic 


pia Anemalistic Lunar Periods are presented in tables. 


PeoultS Suggested by the tidal force plots are discussed 


following 


the tables. 


iompeadeimetanvle the following definitions are used: 


Mag - Magnitudes of earthquakes greater than or 


Se aro tile Magmnmude weomuenedsimethe Cable, 


No - Number of earthquakes in the sample. 


CMe owiie hil-Ssquare Statrstic wath 19 degrees 


of freedom. 


a ZeeCmonc Ct Golme stenitacance level for the 


Vie = 


Ghie-oquare Statistie [Ret Jae 


ites Vi statistic. 


Vn/N - The Vn statistic for Fame wvwaliess on No. 


b = 


ACC = 


Rey 


Polar 


The one per cent significance level for the Vn 
ect 1S 1G. 

NEGepeancemOnntle Milli motihes?s at Ehe one 
Pemeeone Sin teanee slcvc lL. 

Reject On OLericsnlhianypotnesiseat the one 
XSI Gee el cigal se ersyouete | Ike iene 


DijeswOumENc NOGCGUITCNGemOumearthgquakes during the 


Dmdrencmlods are presented an Fagures 3, 4, 5, 6, 7, 8, 9, 


Oe. and» 11. 
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eee CIDAL FORCE COMPONENTS WITH RESPECT TO THE 
SAN FERNANDO VALLEY DATA 


Pugeime IZ 1S an illustration of the tidal force plots 
Heed in the analysis. The occurrence of the earthquake 
in the plot was May 19, 1971 at 0451 GMT. The magnitude 
of the earthquake was 2.6. The radial component of the 
tidal force, as measured in microgals, was positive and 
decreasing in value (negative slope). The north component 
Was negative and increasing in value (positive slope). 

The east component was negative and decreasing in value 
(negative slope). 

The same 110 SFVD earthquakes which were analyzed with 
respect to the Synodic Lunar Period were investigated as 
miesevent presented above. The results are tabulated in 
Table 10. The time frame in which these earthquakes 
meewnrea Was trom April 25, 1971; to October 19, 1971. 

Ten thousand earthquakes uniformly distributed in time 
mere Simulated during the previously mentioned time frame 
and in the same area as the SFVD. These simulated earth- 
Mtekes were analyzed in the same manner. The results are 
also tabulated in Table 10. 

Each component at the time of the earthquake was either 
hecitave or negative in value. The derivatives of each 
of the three components was either positive or negative 
Miemene tame Of the earthquake. Sixty-four combinations of 
the three components and their slopes were then considered. 


inetaple 10 the following definitions are used. 
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Mitre tadial component of th@ tidal force, 


[Oroeomelic NOrEN Component of the tidal force. 
Exowece) Ine Cast component of the tidal force. 
REC - The number of earthquakes recorded with the 


Same combination. 

PREC - The percentage of the earthquakes recorded 
with the same combination. 

SIM - The number of simulated earthquakes recorded 
With the same combination. 

PolM - The percentage of the simulated earthquakes 
with aie Same combination. 

+ “component of tidal force was positive, 

- PECoMmpoment of tidal force was negative. 

P Zvaliie or the component was increasing. 


n - Value of the component was decreasing. 


Miemdaca in Table 10 were used to compare the distribu- 
meen ot Che earthquakes with the distribution of the simulated 
fPerthnaquakes by means of a Contingency Table [Ref. 10]. 

The null hypothesis was that the two distributions are the 
same. The combinations of the three components of the | 
mudaletarce vector and the three derivatives, for which 
mieewaectial and the simulated earthquakes never occurred, 
Were Climinated. This resulted in a 2X34 Contingency Table. 


iicmeir—-oquare Statistic with 35 degrees of freedom was 


comoamted (39.972) and the null hypothesis was accepted. 
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Wee locCUusoolON OF RESULTS AND CONCLUSIONS 


The occurrence of earthquakes, at the larger sample 
mizes. was rejected as a uniform distribution with respect 
to the lunar periods, in all cases except one (Table 8). 

On certain days within the time frames of SFBAD-I and 
SFBAD-II, a larger number of earthquakes occurred relative 
to other days in these time periods. This "swarming" 
effect of numerous earthquakes occurring within hours and 
Mces Of Cach other may indicate the lack of independence 
between the events. The ''swarming" effect could also indi- 
cate that earthquakes do occur in a non-uniform manner. In 
Order to decrease the dependence between the events, if 

the dependence does exist, the data were analyzed according 
wo the Magnitudes of the earthquakes. This method elin- 
Peet cdtne Swarming’ effect, An example of the "swarming" 
effect was on June 10, 1969, when 15 earthquakes occurred. 
When magnitudes greater than or equal to 2.0 were examined, 
Picmniuniper ior Carthquakes which occurred on June 10, 1969, 
Mai-mnecaduceds to five. <The null hypothesis was still rejected 
Pmenomenrce liar periods. The fact remains that the 
null hypothesis was rejected in all cases but one at the 
fancermesamule S1IZes. 

[Meterelabereseiic to Cconpare the polar plots of the 
Synodic Lunar Periods of the SFVD and the SFBAD-II. The 
Umeeostmpetreentage Of Canrtnquakes occurred, in both cases, 


in the first interval, or at the time of the new moon. 
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Also, the time periods of the two sets of data were similar. 
Hewever, thespolareplot of the SFBAD-J was not similar, but 
Meither was the time period. 

imie occttrrence of carthquakes with respect to the 
Praconic™unar Period in the SFBAD-I and SFBAD-II are not 
Biemsanme, When analyzing the SFBAD-I, the null hypothesis 
meomrejected at all magnitudes, but in analyzing the 
SFBAD-II the null hypothesis was accepted at all magnitudes. 

The occurrence of earthquakes with respect to the Anoma- 
ioaere Month in the SFBAD-] and SFBAD-II are similar to 
eachother. Ihe rejection of the null hypothesis at the 
SiiemPerecciu Significance level was the same for the larger 
aanple Sizes. 

The null hypothesis was rejected when analyzing the 
Sev with respect to the Synodic Lunar Period at the one 
per cent significance level when using the Chi-Square Test 
but was accepted using the Vn Test (Both tests would have 
megreercamenc null hypothesis at the 2.5 per cent signifi- 
@mecelever), Ihe rejection was due mainly to the occur- 
emcee: i / Of the 110 earthquakes near the time of the 
new moon. In comparing these 17 earthquakes with respect 
to the combinations of the three components of the tidal 
Pomeemlasmatscissed in section [1ll), nothing of interest 
ciismeemme revealed, Four of the 17 earthquakes had the 
Same combination but the four also occurred on the same 


day. The remaining 15 earthquakes had different combinations. 
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iiemaccepedmee 1m occt1on T1l; that the distributions 
of the 110 actual earthquakes and the 10,000 simulated 
earthquakes were the same seemed quite significant. Further 
analysis into why earthquakes never occur with certain 
combinations and why large numbers of earthquakes occur with 
Other combinations may prove or disprove that the tidal 


momee Cal trigger earthquakes. 
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V. THE COMPUTER PROGRAM 


liicew time )(GMij@e: eaeh origim of ~Ghe lunar periods 
was read into the computer by the year (NA), the month (NB), 
the day (NC), the hour (ND), the minute (NE), and the 
me@omd (St). As discussed in Section JI, the time of each 
euaciaewas converted [Ref. 14] to the corresponding Julian 
Mate (ORGP(1I,J)). 

The time (GMT) of each earthquake was read into the 
computer by the year (MA), the month (MB), the day (MC), 
eacmmeur (MI), the minute (ME), and the second (SE). As 
discussed in Section II, the time of each earthquake was 
converted [Ref. 14] to the corresponding Julian Date 
Pvcwrl)), subsequently each of the JDGMT(1) was con- 
me tcdscomene residual (RESU(1,J)) for each lunar period. 

ice uoliewing list defines the terms used in the Lunar 
Per10d Computer Program: 

NBOX [ies iunben Ot colin tervalseon the 


MES well, (ORI) 


ORGP (1,1) The Voracin seme che ee Synodic Month. 

ORGP (2,1) The Son: Mostiee : cele qe Draconic Month. 
ORGS | ) Iie. Or roonmion s tine (ee Anomalistic Month. 
Peel, 1) INE fore aueye! “oe ielave ie Synodic Month. 
PBR. |) The period of the pele Draconic Monten. 
EeivGs. | } ies werrod O41) ihe 1 Anomalistic Month, 
JDAY Subroutine to convert year, month and day 


to the Julian Day Number [Ref. 14]. 
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BEG (1) 


BEG (2) 


PEGS) 


EN (1) 


EN (2) 


EN(3) 


JDGMT 
2 DLP 


RES 


NCNT 


NPTS 


EXVAL 


eG 


CHISQ 


SHELL 


Mew atrSteorilocim of the Synodic Month 
tO OCCUr Within the time frame of the data. 
(Mice arsenoml cine. thes Draconic Month to 
Ceci witiinmeticmurme <rane ot the data. 
The first origin of the Anomalistic Month 
POMOCCIiM i Miintnic mene Grane oc the data. 
ine las teonreinmear the —Synodic Month to 
Scott ae liimbilcomElnc mcrae Ot etme data, 
hie, lasSteoricinpenrerhce Draconie Month to 
(ecole ELiNc mw hidine Of the data . 
ie las Oto timOonmer nc inonmalistic Month 
POMOCGUDEV itp i-ncmiamne hrane Of the data. 
The Julian Date each earthquake occurred. 
(ie herpes wiiirehince Tenmeanttngqudke occurred 
within the period. 

icc orata ec Oonplincdutor scaGh car tigiuake 
and placed on the unit circle. 

The number of occurrences within each 
Sud) ee male 

The total number of earthquakes. 

Mice te ercomNlulnotmOnmmcaringquakes wi than 
each subinterval. 

They perecentaste Oz wene nuUMbeT Of Occurrences 
Wel il) SeVelon Ssplsnecre sh] a 

Thee Chi-Soquare statistic. 
SUpPOnemicrtOmocOnemenc residuals in 


asecencdmine vOlde Tamineragn 2 |. 
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DIFF Thesvn statistic. 


ROOTS The Vn statistic times YN, 
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